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It is increasingly important to understand, in detail, two-pion correlations measured in p + p and d +A col- 
lisions. In particular, one wishes to understand the femtoscopic correlations, in order to compare to similar 
measurements in heavy ion collisions. However, in the low-multiplicity final states of these systems, global 
conservation laws generate significant A'-body correlations which project onto the two-pion space in non-trivial 
ways and complicate the femtoscopic analysis. We discuss a model-independent formalism to calculate and 
account for these correlations in measurements. 
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I. INTRODUCTION 

The unique and distinguishing feature of heavy ions is their 
large (relative to the confinement scale) size and the possibil- 
ity to generate bulk systems which may be described in ther- 
modynamic terms, allowing to discuss the Equation of State 
of strongly-interacting matter. The primary evidence for the 
creation of bulk matter at the highest energies is the existence 
of strong collective flow [1]. The dominant feature of flow is 
the correlations between space and momentum which it gen- 
erates; thus, momentum-only observables such as pj spectra 
and azimuthal anisotropies |l2|, |3t |4j, |5|] represent only an indi- 
rect projection of the effect. Femtoscopic measurements ac- 
cess space as a function of particle momentum, thus providing 
the most direct probe of the most crucial feature of heavy ion 
collisions [c.f. e.g. 6J. In particular, flow is manifest by a neg- 
ative correlation between the "HBT radius" and the transverse 
mass (mj) of the particles 101. 

Clearly, then, a detailed understanding of femtoscopic mea- 
surements in heavy ion collisions is crucial to proving the ex- 
istence of, or probing the nature of, the bulk system gener- 
ated in the collision. It is in fact possible to quantitatively 
interpret both the femtoscopic and momentum-only observa- 
tions at RHIC in consistent, flow-dominated models of the 
system [e.g. 8]. All seems well. 

However, two-pion femtoscopic measurements are also 
common in e+ + or p + p{p) collisions |9]. In these 
collisions, too, "HBT radii" are observed to fall with mj. 
Speculations of the physics behind this observation have 
included Heisenberg uncertainty-based arguments, string- 
breaking phenomena, and temperature gradients; an excellent 
overview may be found in [ 10]. Typically, however, one might 
not expect the system created in a p + p collision to exhibit 
bulk behaviour similar to that from heavy ion collisions. 

Quantitative comparisons between femtoscopic measure- 
ments in A + A and p+ p systems have been complicated be- 
cause techniques for event-mixing, frame definitions, and the 
like, have been different in the particle-physics and heavy-ion 
communitites. As importantly, kinematic acceptance and col- 
lision energies are usually quite different. Recently, however, 
the STAR experiment has reported the first preliminary study 
of directly-comparable femtoscopic measurements from A + A 
and p + p systems llTl at the same ^/snn, using the same de- 
tector, and with identical techniques. The results indicate that 
the femtoscopic probe of flow- falling "HBT radii" with mj- 



is essentially identical in the small and large systems. This 
might signal an unexpected "universality" in the spatial sub- 
structure of hadronic and heavy ion collisions. Unravelling 
the physics behind this similarity might provide new insight 
into the dynamical space-time substructure of both hadronic 
and heavy ion collisions. 

Before drawing strong physics conclusions from "HBT 
radii" coming from fits to the pion correlations measured in 
p+ p collisions, however, the measured correlation functions 
themselves must be understood in detail. The STAR data 
show clear non-femtoscopic correlations which must be dis- 
entangled from the femtoscopic ones flT]. Femtoscopic cor- 
relations are those which depend directly on the two-particle 
separation distribution [c.f.|6[]. Non-femtoscopic correlations 
may arise from string or jet fragmentation, resonance decay, 
or global conservation laws. 

In this work, we explore the projection of A^-body Energy 
and Momentum Conservation-Induced Correlations (EM- 
CICs) onto a two-particle relative momentum correlation 
function. In Section we briefly discuss the harmonic rep- 
resentation of the correlations which best illustrates the effect. 
In Section|III]we discuss EMCICs generated by a Monte Carlo 
event generator containing only global conservation laws. A 
method to calculate analytically (but using distributions from 
the data) EMCICs is shown in Section |IV] This provides an 
"experimentalist's formula," given in Section[Vl useful to dis- 
entangle EMCICs from the data, allowing a femtoscopic anal- 
ysis of the correlation functions. We summarize in SectionrvTl 



II. SPHERICAL HARMONIC DECOMPOSITION OF 
CORRELATION FUNCTIONS 

At asymptotically high relative momentum |^| (or \k*\), 
femtoscopic contributions to the the correlation function 
(those described by the Koonin-Pratt equation [discussed in 
6]) must approach a constant value, usually normalized to 
unity, independent of the direction of q. Preliminary STAR 
measurements lITlll of small systems. Figure [T] show clear 
non-femtoscopic correlations in addition. Also shown is a 
fit with the commonly Gaussian (with Coulomb suppression) 
functional form |@]. Clearly, the fit is a poor representation of 
the data. We stress, however, that it is not the (non-)Gaussian 
nature of the source at issue here; any source function will 
lead to vanishing femtoscopic correlations at large |^| and will 
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FIG. 2: Preliminary STAR two-pion correlation functions fTHl pre- 
sented in the SHD representation. 



FIG. 1 : Preliminary STAR two-pion correlation functions ITTIl pre- 
sented as ID projections in the Bertsch-Pratt decomposition. 



thus contradict the data. 

We further stress that the problem is not one of normaliza- 
tion. Shown in the Figure is the common representation of the 
3 -dimensional correlation function into three 1 -dimensional 
axes [cf The projections, then, are not independent and 
cannot be independently normalized. The problem is that the 
value approached at large \q\ depends on the direction in q 
space. 

One-dimensional projections present a limited tool for ex- 
ploring detailed structure of the three-dimensional correlation 
function. The spherical harmonic decomposition (SHD) is a 
much more efficient representation of the data which uses all 
of the data to show the shape of the correlation function in 
3D q space. There, the spherical harmonic coefficients A/ 
which depend on 2 = |^|, are calculated as 



A/,„,(e)= c(e,cose,-,^,)- 

bins / 

Y,,,n (cOSe;,(|),)/7.,„ (cOS0/,Acos9,A.|,) , (1) 

where F/ ^ represents a numerical factor correcting for finite 
bin sizes AcosO ™d ^if\ it turns out not to depend on (|),-. The 
angles 9 and (|) are related to the Bertsch-Pratt Cartesian coor- 
dinate system through 

?o = 2sin9cos(|), ^.s = 2sin0sin(|), qi^Qcos,Q. (2) 

See lfl2ll for a complete discussion. 

Preliminary STAR correlation functions in the SHD repre- 
sentation fll'l are shown in Figure|2] Coefficients for / > 4 are 
much less significant, compared to errorbars; to good approx- 
imation, the non-femtoscopic behaviour is quadrupole (I = 2) 
in nature. 

The presence of non-femtoscopic correlations is clear from 
the non-vanishing behaviour of A/^o,m's at large Q. However, 



it is by no means clear that these contributions to the cor- 
relation function are confined to large Q. Thus, one cannot 
attempt to interpret the low-Q region only in terms of fem- 
toscopic coiTelations, while parameterizing or ignoring the 
large-Q region; see lll3[l for further discussion. 



III. EMCICS GENERATED BY THE GENB OD MONTE 
CARLO GENERATOR 

Non-femtoscopic correlations may arise from a variety of 
sources. Jets will clearly induce momentum-space correla- 
tions between its fragmentation products. While this cannot 
be discounted, the low momentum of the pions under consid- 
eration ipT ^ 0.2 GeV) puts us squarely in the region in which 
factorization breaks down and the jet interpretation becomes 
significantly murkier We do not explore this possibility here. 
In the kinematic region under consideration, string fragmen- 
tation may play a role; this is an area for future study, though 
significant model-dependence will be present. 

Resonances induce correlations among daughters as well; 
while these might even dominate 7i+7i^ correlations, they 
should be negligible for identical pion correlations. Collec- 
tive bulk flow (e.g. anisotropic elliptic flow) will generate A^- 
body correlations which will project onto the two-body space. 
Non-femtoscopic correlations of the type observed by STAR 
in small systems are not, however, observed in AuH-Au colli- 
sions, despite the fact that elliptic flow is much larger there; 
therefore, we do not believe that collective flow generates the 
observed effects. 

Without doubt, one physical effect which must be at play 
is momentum and energy conservation. As global conserva- 
tion laws, these provide an A^-body constraint on the event, 
which projects down onto 2-body spaces. The observed non- 
femtoscopic effects il ill become more and more significant as 
the multiplicity (A^) of the event decreases, as expected from 
conservation laws. It is these EMCICs which we focus on 
here. 
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phase-space weight = 5.62509741 e-36 




FIG. 3: A high-probability muhiphcity-30 event calculated by GEN- 
BOD . Lines correspond to particle momenta Px,Py,Pz- 



phase-space weight = 3.78350585e-44 




FIG. 4: A low-probability multiplicity-30 event calculated by GEN- 
BOD . Lines correspond to particle momenta Px,Py,Pz- 



To clearly understand the role of EMCICs, we would like 
to have events in which there is no other physics involved be- 
sides the conservation laws. Such a tool has been provided 
almost 40 years^ ago in the form of the GENBOD computer 
program [see 14, for an excellent writeup of the method and 
physics] in the cern library. Given a requested total multi- 
plicity (A^), a list of masses (m,) of emitted particles, and a 
total amount of energy (fitot) to distribute among them, GEN- 
BOD returns an event of random momenta (four vectors pj), 
subject only to the condition of energy and momentum conser- 
vation. More importantly, it returns, for each event, a weight 
proportional to the probability that the event will actually oc- 
cur in nature. This weight is proportional to the phasespace 
integral Rn 



R^ 




(3) 



where P = i^uOj is the total momentum four- vector of the 

event. See fl4j for a practical iterative prescription to calcu- 
late Rn. Thus, it is a much different tool than, say RQMD, in 
which each event returned may be treated as equally probable. 

We select (via monte-carlo) GENBOD events according to 
their weight and run them through identical software as used 
for experimental analysis. Fortunately, the code is fast, since 
one must calculate large statistics from which to select. This 
is because the phase-space weights vary by large factors. As a 
very extreme case, Figures[3]and|4]show a likely and unlikely 
event, respectively, for multiplicity A^ — 30. As one would 
expect, the "rounder" event is more likely, though one might 
be surprised by the factor of a hundred million between the 
probabilities. 

Figures |5] and |6] show the A/ „,'s calculated by GENBOD for 
18-pion events without and with a selection of |r|| < 0.5, re- 
spectively. Note that this cut applies to the pions which are 
used in the analysis, not to the set of particles for which en- 
ergy and momentum is conserved; energy and momentum is 
always conserved for the full event. Clearly visible are signif- 
icant and nontrivial A/„,'s due only to EMCICs. We observe 
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FIG. 5: SHD coefficients for GENBOD -generated events consisting 
of 18 pions, as measured in the pair CMS frame. Green squares are 
A/ „,'s from the GENBOD events. For discussion of the other symbols, 
see SectionlTVl 



also that the / = 4 coefficients are about an order of magnitude 
smaller than the I = 2 ones; this is generically expected [cf 
12]. Comparing the two figures, it is clear that kinematical se- 
lection has significant effect on the EMCIC effects. Also sig- 
nificant (but not shown) is whether one includes other species 
(say protons) into the mix of emitted particles. 

Comparison of Figures|7][8]and|9]makes clear the multiplic- 
ity dependence of the EMCICs. As expected, lower multiplic- 
ity events show a greater effect. Also (not shown), increasing 
the amount of energy to be distributed among the particles, for 
fixed multiplicity, decreases EMCICs, as one expects. 

Finally, we note that EMCICs can affect the correlation 
function even down to very low Q, again reminding us that 
we cannot (responsibly) ignore these effects in a femtoscopic 
analysis. 
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FIG. 6: Same as in Figure|5] except only using pions with |r|| < 0.5 
in the correlation function. 
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FIG. 8: Same as Figure|7] but for 9-pion events. 
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FIG. 7: SHD coefficients for GENBOD -generated events consisting 
of 18 pions, as measured in the pair LCMS frame. Green squares are 
A/ „/s from the GENBOD events. For discussion of the other symbols, 
see SectionlTVl 



IV. ANALYTIC CALCULATION OF EMCICS 

Now then, EMCIC effects generated by GENBOD "resem- 
ble" the experimental data, but it is likely unwise to use GEN- 
BOD itself to correct the data for several reasons. Firstly, there 
is strong sensitivity to the (not completely measured) number 
and species-mix of all particles emitted in the event, includ- 
ing neutrinos and possible magnetic monopoles (or, less exoti- 
cally, particles escaping detector acceptance). Secondly, there 
is strong sensitivity to the energy "available" in the event; it 
is not obvious that this is ^/snn of the collision. Clearly, EM- 
CIC effects depend on the individual momenta pi and p2 of 
the particles entering the correlation function. This will de- 
pend on acceptance, efficiency, kinematic cuts (both purpose- 
ful and those imposed, e.g., by particle-identification limita- 
tions), and, to a degree, the underlying single-particle phas- 
espace. (While correlation functions are insensitive to the 
single-particle phasespace, the correlations which they mea- 
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FIG. 9: Same as Figure[7] but for 6-pion events. 



sure may, in fact, depend on this phasespace, due to physical 
effects.) 

Thus, one would like to calculate EMCICs, based on the 
data itself. An excellent start in this direction has been pre- 
sented by Danielewicz and later by Borghini, Dinh and 
Ollitraut 1 16], in which they considered the effect of EMCICs 
on two-particle azimuthal correlations (elliptic flow V2). They 
considered transverse momentum (Pj) conservation only, but 
the extension of their formula to three, and even four lUSll di- 
mensions is straightforward. 

We consider the single-particle distribution of a mass-m, 
particle Mnaffected by EMCICs 



f{p,)=2Ei — 



(4) 



The ^-particle distribution (k less than total multiplicity A^), 
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including EMCICs is then 



fc{p\,...,Pk) 



n/>') 



/ {y\%,+, (p] - m^) / {pi)) 5^ (It, pt - P) 

J [UU d'Pj?> {Pj - ^) f iPj)) §^ (It 1 Pi - P) 



(5) 



(Note the difference between numerator and denominator in 
the starting value of the index j on the product.) 

According to Equation |5] the fe-body momentum distribu- 
tion, including EMCICs, is the A;-body distribution not af- 
fected by EMCICs- i.e. just an uncorrected product of 
single-particle distributions- multiplied by a "correction fac- 
tor" which enforces the EMCIC. The numerator of this factor 
just demands that the remaining N — k on-shell particles are 
configured so as to conserve total energy and momentum, and 
the denominator just normalizes the distribution. 

Using the central-limit theorem (valid for "large enough" 

— k), we find 

fc{pi,-,Pk) = (jJfiPi)] X 



N 



N-k 



■ exp 



^ {LLip..,-{p,))) 

2{N-k)al 



where 



and 



ol^{pl)-{p,)' 



(pD^ 



dpf{p) ■ Pu- 



(6) 



(7) 



(8) 



Naturally, (/J(;j= 1,2,3)) — 0- (In these equations, we now as- 
sume only one species of particles, so that no species label 
is needed for (p^). This is only for simplicity of notation 
here; results, including the "experimentalist's formula" below, 
only become more cumbersome to write, but are similar oth- 
erwise.) 

Note that even the single-particle momentum distribution is 
affected by EMCICs 

fc{Pi)=f{pi)-i^\ X (9) 



exp 



1 



N- 1 
Ply 



9?^ (Ei-{E)f 



The fc-particle correlation function is defined as the mea- 
sured (i.e. EMCIC-affected) ^-particle yield divided by the 
product of the measured single-particle yields 

_ fc{ph-;Pk) _ {wmf ^, 



C{pi,...,pk) ■ 



fc{pi)---fc{pk) {^f 



(10) 



exp 


2{N 


-k) {^f=l ( 




, {l.\{E,-(E)))-\ 


exp 


-1 s^k 

2(N-l) 




(E,-(E)f V 





An important point: EMCICs result from the constraint that 
the event's energy-momentum is the same fixed number for 
all pairs in the event. This is true in the laboratory frame, but 
not in LCMS or pair rest frame. Thus, while one may bin the 
correlation function in the frame of one's choice, the momenta 
which appear on the right-hand-side of Equation [TOl must be 
in the laboratory system. 

To first order in the two-particle correlation function 
becomes 



1 



C{p\,P2) = 1 - 
,Pl,T-P2.T , P\,z-P2,Z {El 



(11) 



{E)){E2-{E)) 



(pD 



{E^)-{Ey 



where we have taken (p^) = (p^) = (pj) /2. In what follows, 
we shall refer to the first, second, and third terms within the 
parentheses of Equation [TTI as the "pr term," "/?,. term," and 
"£■ term," respectively. 

If we know A^, (/j^), (p^), {E^), and (E) from the data, 
we can calculate EMCICs using Equation[TO] Better yet, if 
is large enough, then we can use Equation [TT] This is what 
is done in Figures |5]|8] The black circles, blue stars, and red 
triangles show the pj, Pz and E terms, respectively, from the 
first-order expansion (EquationfTTTi. while the open circles and 
orange inverted triangles represent the results of Equation [TTI 
and Equation [Tol respectively. 

Several observations are in order. Firstly, each of the three 
terms in Equation [TT] produce non-trivial behaviour of the 
A/„,'s, interfering with each other in interesting ways. We find 
also that the p- term affects A2,2\ this was initially surprising 
since A2,2 quantifies the behaviour of the correlation function 
in the "out-side" plane, while z is the "long" direction in the 
Bertsch-Pratt system. 1 19] Clearly, EMCICs projected onto a 
2-particle space are non-trivial objects. 

It is seen that the first-order expansion (Equation[TTli agrees 
well with the full expression (Equation[TOli well for A^ >~ 10. 
Such multiplicities are relevant for the p + p measurements 
done at RHIC (especially recalling that A^ includes all parti- 
cles, even unmeasured ones). We see also that the analytic cal- 
culations (open circles and inverted triangles) approximate the 
results of the GENBOD simulation (green squares), especially 
as the multiplicity and total energy of the event increases; in- 
creasing agreement for large A^ and is expected, given the 
approximations leading to our analytic expressions. We ob- 
serve also that the analytically-calculated expressions respond 
identically to the kinematic cuts as does the simulation (c.f. 
Figures [5] and [6]l. 

Never do the analytic calculations never reproduce exactly 
the simulations; we discuss this further in the next Section. 



V. AN EXPERIMENTALIST'S FORMULA 

Even for large A^ and energy, the calculations do not ex- 
actly reproduce the EMCIC effects in the simulation. One 
reason for this may be found, in fact, in the definition of the 
average values (e.g. {p^}) themselves. In Equation[8] average 
quantities are calculated using the distribution f{p), which 
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is not affected by EMCICs. Naturally, the only measurable 
distribution available to the experimentalist (even when GEN- 
BOD simulations serve as the "experiment") is fc{p)- 

Thus, it appears the experimentalist cannot plug her data 
into the equations |7] [8] and [TT] to fully calculate EMCICs. 
However, such an ambition would have been hopeless any- 
how. After all, even the total multiplicity (again, including 
photons etc) is rarely fully measured. And finite kinematic 
acceptance (e.g. in r|) will require extrapolation to calculate, 
e.g. {p^). 

To the practicing femtoscopist, there is a natural solution. 
Having at hand (1) educated guesses for the quantities A^, {E^) 
etc, and (2) a physically-motivated functional form which 
connects these quantities to the correlations we'd like to un- 
derstand, we perform a fit. Let us rewrite EquationfTTIas 



C{p\,P2) = 1 -Ml ■ {pi,T ■ P2,t} -M2- {pi,z- P2,z} (12) 



^M3 • {El ■E2}+M4 ■ {El +E2}-^. 



where 



Ms 



M, = 



N{pI] 



, M2 



1 



N{{E^)-{E)^) 



- N{pI 
M4 



N{{E^)-{Eyy 



(13) 



The notation {X} in Equation [12] highlights the fact that X 
is a two-particle quantity which depends on pi and p2 (or q, 
etc): {X} [q). From a practical point of view, X is averaged 
over the same q bins as used for the correlation function. This 
involves nothing more than adding four more histograms to 
the several akeady being constructed by the experimentalist 
as she runs through all pairs in the data. The binned functions 
{X{ then automatically reflect the same event and particle se- 
lection as the correlation function. 

It is appropriate here to re-emphasize the point made in ref- 
erence to Equation [TOl The event's total energy and momen- 
tum is a fixed quantity in a fixed (e.g. lab) frame; in partic- 
ular, the momentum in Equation [TOl is assumed P = 0- i.e. 
the collision-center-of-mass (CCM) frame is assumed. In a 
pair-dependent frame (e.g. pair center-of-mass PCM or lon- 
gitudinally co-moving system LCMS), the event's energy and 
momentum will depend on the pair. EMCICs, therefore, must 
be calculated with CCM momentum. Thus, in the function 
{p\,z ■ Pi.z} Pi.z and P2.Z must be calculated in the CCM 
frame, while the binning variable q should be in whatever 
frame one chooses to construct the correlation function C. 

The parameters M, defined in Equation [13] on the other 
hand, are global and independent of pi and p2. It is these 
which we will use as fit parameters. The task is then fast and 
straightforward; the EMCIC part of the correlation function 
C{q) is simply a weighted sum of four functions. Indeed, one 
may calculate coefficients as in Equation [T] for the four new 
functions. For example 



Al%(Q)^ E {/'l,z-/'2,J(fi,C0Se,-,C^,) 



(14) 
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FIG. 10: A/ „/s from 18-pion GENBOD -generated events. Green in- 
verted triangles (often underneath black squares) is the correlation 
function (measured in PRF) from GENBOD . Filled brown circles, 
filled blue stars and filled red triangles show, respectively, the "pj-," 
"p^," and terms, defined in Equation 1 1 1 [ black filled squares 
show their sum. Open symbols of the same shape and color (identi- 
fied as "FIT" in the legend) show corresponding terms, except with 
weights (see Equation I73]l adjusted to maximize agreement between 
the open black squares and the simulation. 



etc. Then, thanks to the linearity of Equation [12] and the 
orthonormality of F/ ,„'s, the measured A/.;„'s themselves are 
similarly just weighted sums of harmonics 

A/,,„(e) - 5,.o • (1 -MI/M3) -Ml -All^ (Q) (15) 



-M2 -A'Jl {Q) -M3 . A{^f ) {Q) +M4 . Ar„r^ [Q) . 

Treating Equation[T5]as a fit, we have a few (say six, for I < 4) 
one-dimensional functions to fit with four adjustable weights. 

A first example of such a fit is shown in Figure [TO] Again 
the GENBOD simulation is compared to the first-order form of 
Equation [12] The filled circles, stars and triangles show the 
"Pt' (Ml), "pz' (M2), and "E" (M3 and M4) terms when the 
weights (Equation [T3Tl are calculated directly from the events, 
as discussed in Section iTVl Treating the M, as adjustable pa- 
rameters leads to a slightly different weighting of the terms, 
and a slightly better fit to the data. 

Lest we forget, our original goal was not to understand EM- 
CICs per se, but to extract the femtoscopic information from 
measured two-particle correlations. Assuming that the only 
non-femtoscopic correlations are EMCICs, one may simply 
add the femtoscopic terms (e.g. Gaussian in {qo,qs,qi) space 
or whatever) to the fitting function in Equation[T2]or[T5] Com- 
mon femtoscopic fitting functions usually contain ~ 5 param- 
eters (e.g. Rout) In the imaging technique iITtIi . one assumes 
the separation distribution is described by a sum of splines 
(rather than, say, a Gaussian); here, too, there are usually 4-5 
fit parameters (spline weights). We have found that the num- 
ber of fit parameters now must be doubled to account also for 
EMCICs. This is a non-trivial increase in analysis complexity. 
However, we keep in mind two points. 
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Firstly, the increased work is neccessary. EMCICs (and 
possibly other important non-femtoscopic correlations) are 
present and increasingly relevant at low multiplicity. One op- 
tion is to ignore them, as has sometimes been done in early 
high-energy experiments. However, with the new high-quality 
data and desire for high-detail understanding at RHIC, ignor- 
ing obvious features such as those seen in Figures [T] and |2] 
is clearly unacceptable. Perhaps a slightly better option is to 
invent an ad-hoc functional form with no real physical basis 
(and often manifestly wrong symmetry [cf jTst]), which intro- 
duces new parameters in any case. We hope that the results 
here present a relatively painless, and considerably more re- 
sponsible, third option. 

Secondly, while the non-femtoscopic EMCICs are not con- 
fined to the large-2 region (an important point!), the femto- 
scopic correlations are confined to the small-Q region. There- 
fore, one hopes that the addition of four new parameters 
to the fit of the correlation function will not render the fit 
overly unwieldy. While we can not expect complete block- 
diagonalization of the fit covariance matrix, one hopes that the 
Mi are determined well enough at high Q that the femtoscopic 
fit parameters can be extracted at low Q. 

VI. SUMMARY 

To truly claim an understanding of the bulk nature of matter 
at RHIC and the LHC, a detailed picture of the dynamically- 
generated geometric substructure of the system created in 
heavy ion collisions is needed. It is believed that this substruc- 
ture, and the matter itself, is dominated by strong collective 
flow. The most direct measure of this flow is a measurement 
of the space-momentum correlation (e.g. R [mT]) it generates. 
The physics of this large system, and the signals it generates, 
should be compared to the physics dominating p + p colli- 
sions, as is increasingly common in high-y^j- studies at RHIC. 
For the small systems, however, non-femtoscopic effects con- 
tribute significantly to the correlation funcion, clouding the 
extraction and interpretation of the femtoscopic ones. 

EMCICs, correlations generated by kinematic conservation 



laws, are surely present and increasingly relevant as the event 
multiplicity is reduced. Using the code GENBOD to study cor- 
relation functions solely driven by EMCICs, we found highly 
non-trivial structures strongly influenced by event characteris- 
tics (multiplicity and energy) and kinematic particle selection. 

We extended the work of Danielewicz and OUitrault to in- 
clude four-momentum conservation and applied it to corre- 
lation functions commonly used in femtoscopy. We found 
structures associated individually with the conservation of 
the four-momentum components, which interfere in nontrivial 
ways. Comparison of the analytic EMCIC calculations with 
the GENBOD simulation gave confidence that the approxima- 
tions (e.g. "large" multiplicity A^) entering into the calculation 
were sufficiently valid, at least for multiplicities considered 
here. We further showed that the full EMCIC calculation can 
safely be replaced with a first-order expansion in 

Based on this first-order expansion, we developed a practi- 
cal, straight-forward "experimentalist's formula" to generate 
histograms from the data which are later used in a generalized 
fit to the measured correlation function, including EMCICs 
and femtoscopic correlations. 

The huge systematics of results and interest in femtoscopy 
in heavy ion collisions is renewing similar interest in the 
space-time signals from p+ p collisions. Direct comparisons 
between the two systems are now possible at RHIC and have 
already produced intriguing (albeit preliminary) results. Very 
soon, p + p collisions will be measured in the LHC exper- 
iments, and the heavy ion experimentalists will be eager to 
apply their tools. The femtoscopic tool is one of the best in 
the box- so long as we keep it sufficiently calibrated with re- 
spect to non-femtoscopic effects increasingly relevant in small 
systems. 
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